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ABSTRACT 

Aims. For the past 25 years, BW Vulpeculae has been the topic of period analyses centered on a secular period change with a periodic 
variation superposed, presumed to be due to light time effects in a binary system. According to this paradigm, one would expect what 
seems like a period increase of about 0.5 s during or soon after the year 2001. 

Methods. I have continued photometric monitoring through the year 2012, adding 35 new timings of maximum and minimum light. 
Results. This expected change in period did not occur, which rules out that interpretation of the period variation. As of 2012, the 
observed timings are about two hours early compared to those predicted by the quadratic ephemeris, but are very close to those 
predicted by the linear ephemeris. 

Conclusions. In fact, the period has remained constant for the last 32 years, indicating that the previous epochs of constant period are 
almost certainly the correct interpretation, though the cause of the period changes is still not clear. Continued photometric monitoring 
of BW Vul leads to the conclusion that the period changes are abrupt, followed by epochs of constant period lasting between 12 and 
at least 32 years. 
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1. Introduction 

Much has been written in the last 30 years about changes in 
the period of BW Vulpeculae, the monoperiodic /3 Cep star with 
the largest known amplitude (van der Linden & Sterken [19871 ). 
Petrie d 19541 first suggested a constant rate of period increase 
dP/dt = +3.7 s/cent based on spectroscopy, while Cherewick & 
Young ( 1975 ) confirmed this with photometry but derived a rate 
half as large. These large positive dP/dt would indicate that, if 
due to the evolution of the star, BW Vul must be in the shell 
hydrogen burning phase. This contradicts some evidence that /3 
Cep stars in clusters are in the late core burning stage, where the 
rate of period change is well less than 1 s/cent, and it certainly 
rules out the overall contraction phase where the period should 
decrease. 

Odell (1984) noticed what appeared to be a periodic varia- 
tion superposed on the quadratic ephemeris with half-amplitude 
about 17 minutes and about 25 year period, and suggested ei- 
ther a light-travel time effect (LiTE) of a small-mass companion 
or two pulsation modes beating with that long period. Pigulski 
( 11993b solved for the binary orbit postulated by Odell, and found 
a mass function of 0.012 Mq, with a period of 33 years and an 
eccentricity of 0.46. For reasonable assumptions about the incli- 
nation and mass of the primary, this suggests the secondary is 
less than 2.5 Mq, so it is undetectable in extant observations. 

Tunca (119781) proposed that the period of BW Vul had been 
constant, but abruptly increased in 1972 by about 0.5 s. Similar 
interpretation was made by Chapellier (119851 ) who suggested 
abrupt period changes in 1931 of +0.61 s, 1945 of +0.50 s, and 
1971 of +0.52 s, with the period changes taking place over sev- 
eral years. He finds that the standard deviation of the residu- 
als for the linear ephemeris fits are reduced by a factor of four 
over that of the single parabolic ephemeris, but this does not take 
into account the putative variation due to LiTE, which Pigulski 
( 119931 shows also substantially reduces the residuals. 



Chape llier's ( 1985 ) suggestion of abrupt period changes fol- 
lowed by epochs of constant period just happens to mimic a con- 
stant dP/dt with a periodic variation superposed on it, but the 
latter makes a very strong prediction that another seeming pe- 
riod increase should have occurred in or around 1980-81, and 
this exactly happened (Chapellier & Garrido 11990b . According 
to Chapellier's ( 1985 ) suggestion, this would be fortuitous coin- 
cidence, though he proposes no physical mechanism for causing 
the period to change, but relies on a possibility that some con- 
vective process in the star is responsible. It is unfortunate that 
the international campaign (Sterken et al. [1986 ) to monitor B W 
Vul was organized for the 1982 observing season, and only one 
timing was done in the year of the period change, 1981. This led 
Chapellier & Garrido (1990) to suggest that both the amplitude 
and timings became unstable for three years during the change. 

All this has caused our understanding of BW Vul to be in 
upheaval, which has been mostly ignored because of the seeming 
predictive power of LiTE, and so the case has been considered 
settled (see eg Horvath et al. |1998l hereafter HGF). BW Vul went 
through more than two complete cycles of LiTE between 1934 
and 2001. Two excellent review papers on the (O-C) (residuals) 
diagram, Sterken ( 120051 ) & Zhou ( [1999] ) both use BW Vul as the 
illustrative example of a star showing LiTE effects. 

However, the LiTE model makes another prediction, that 
there should again be a period change around 2002, with an in- 
crease of about 0.5 s. The purpose of this paper is to report that 
such a change did not happen, and therefore both the large rate 
of change of period and the LiTE can be ruled out, at least in this 
star. The failure of the star to change period again is manifest in 
a discrepancy by 2012 of over two hours compared to the predic- 
tion of the cannonical model, as shown in Section 2. Though this 
does not confirm the piecewise-linear ephemeris with abrupt pe- 
riod changes, it is consistent with it; continued monitoring might 
clarify the situation in this regard. 
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2. Observations 

BW Vul was monitored at Lowell Observatory from 1995 to 
201 1 using three different instruments. Before 2002, the white 
photoelectric photometer was used on the 31 -inch telescope. On 
JD 52948 and 55844, the NASACam CCD was used in robotic 
mode on the 31 -inch telescope (note all HJDs have 2,400,000 
subtracted). The rest of the measurements reported here were ob- 
tained with the photoelectric photometer permanently mounted 
on the 21 -inch telescope. Primarily the Stromgren b and y fil- 
ters were used, though others were employed occasionally. In 
all cases the comparison star was HD198820=HR7996, which 
is similar color and brightness to the variable, so no color terms 
were included. The comparison was used to find the extinction 
coefficient, and all measurements were corrected for extinction; 
the comparison was interpolated to the time of the variable mea- 
sure, and the magnitudes were subtracted. 

These observations are summarized in Table [T] Col. 1 gives 
the cycle number based on the maximum of Huffer (1938) be- 
ing cycle 0. Rather than convert times of minimum, cycles with 
fraction 0.443 were assigned (Sterken et al 1987J- Col. 2 gives 
the HJD of the timing - this was found by folding back in time 
the light curve so the rising and falling branches coincided; the 
time at the fold was taken as time of maximum or minimum, 
thus making obvious any effects of the stillstand. The fourth col- 
umn is the residual, or observed time minus the computed time 
(O-C) based on a linear ephemeris fit to timings between 1982 
and 1998 (Eq. 4); these are plotted in Fig. 4. The fifth column is 
(O-C) based on an ephemeris which includes a quadratic fit (Eq. 
2) between 1934 and 2001 (two complete cycles of the supposed 
LiTE) and a removal of the LiTE (Eq. 3); these are plotted in 
Fig. 3. 

Note that neither ephemeris utilizes the potentially contro- 
versial timings after 2002. The timings used for the fits and the 
graphs include supplemental times from Sterken (priv. comm.), 
HGF, and the AAVSO. The entire list of 335 timings is available 
from the author. 



Table 1. Timings for BW Vulpeculae 



Cycle 


HJD (2400000+) 


Max/Min 


(O-C) 1 


(o-cy 


105517.443 


50015.6330 


Min 


2.2 


3.1 


105523. 


50016.7540 


Max 


8.8 


9.7 


107209. 


50355.7087 


Max 


1.1 


2.9 


108568. 


50628.9300 


Max 


4.9 


6.8 


108572. 


50629.7340 


Max 


4.6 


6.5 


109134. 


50742.7213 


Max 


5.5 


7.3 


111069. 


51131.7371 


Max 


-0.5 


-0.2 


116460.443 


52215.6610 


Min 


8.5 


-3.7 


116461. 


52215.7757 


Max 


13.6 


1.4 


116465. 


52216.5737 


Max 


4.7 


-7.5 


116480. 


52219.5990 


Max 


18.6 


6.3 


116485. 


52220.6055 


Max 


20.4 


8.1 


117362. 


52396.9143 


Max 


10.7 


-4.7 


117615. 


52447.7815 


Max 


15.2 


-1.3 


117616. 


52447.9690 


Max 


-4.3 


-20.8 


117620. 


52448.7823 


Max 


8.8 


-7.6 


117695. 


52463.8650 


Max 


15.2 


-1.6 


120106. 


52948.5758 


Max 


6.4 


-20.9 


120106.443 


52948.6670 


Min 


8.3 


-19.0 


120107. 


52948.7879 


Max 


22.3 


-5.0 


127423. 


54419.6259 


Max 


23.3 


-47.7 


134084. 


55758.7718 


Max 


12.3 


-112.3 


134084.443 


55758.8585 


Min 


7.7 


-116.8 


134188. 


55779.6750 


Max 


4.6 


-120.9 


134188.443 


55779.7662 


Min 


6.5 


-119.0 


134193. 


55780.6825 


Max 


7.8 


-117.7 


134193.443 


55780.7708 


Min 


5.6 


-120.0 


134263. 


55794.7563 


Max 


8.9 


-117.3 


134446.443 


55831.6397 


Min 


12.5 


-115.3 


134447. 


55831.7505 


Max 


11.9 


-115.9 


134511. 


55844.6161 


Max 


10.2 


-118.2 


134511.443 


55844.7042 


Min 


7.7 


-120.8 


134512. 


55844.8123 


Max 


3.2 


-125.2 


135542. 


56051.8905 


Max 


7.5 


-130.4 


135542.447 


56051.9803 


Min 


7.4 


-130.5 



(O-C) 1 is the residual to a linear fit of data between 1981 and 2002. 
(O-C) 2 is the residual to a quadratic fit plus LiTE model (HGF Eq.8) 



3. Analysis 

The method of analysis adopted here is to utilize the exact steps 
of HGF. Theirs is the most recent attempt to understand the pe- 
riod behavior of BW Vul, including data through 1997, and it re- 
produces the results of several earlier investigations. That paper 
also gives all the fitting functions which result from the secular 
period increase and LiTE. 

The first step in analyzing the period of a variable star is to fit 
a linear ephemeris to the timings vs cycle number; the slope of 
this line would be the period. Then, the (O-C) diagram (resid- 
uals plotted vs HJD) will show how well that ephemeris works, 
and whether there is some trend away from it. When this is done 
for B W Vul, the resulting linear equation (HGF Eq. 2, taken from 
Sterken (fT99"3l )[]is: 



HJD miX = 28802.5487 + 0.201038 E 



(1) 



1 I have changed the initial epoch to reflect that I use timings of max- 
imum light for all their equations. This is because most historical tim- 
ings are time of maximum light. I make a plea here to modify the cycle 
count rather than the actual timings, as I've done in Table[T] This, I be- 
lieve, is more academically honest, as it doesn't require modifying the 
data. Likewise, I use actual times of radial velocity zero crossing, with 
a phase offset by 0.522 in the cycle count rather than infer a time of 
maximum or minimum light from RV curves. 




20000 25000 30000 35000 40000 45000 50000 55000 60000 

Heliocentric Julian Date (2400000+) 

Fig. 1. Residuals for HGF Eq 2 (my Eq 1), the solid lines are fits 
of a parabola and five piecewise-linear fits to subsets of the data. 



which gives the residuals shown in Fig. [T] It is clear from that 
figure that a single constant period doesn't represent the timings 
over the entire dataset (there are residuals up to +400minutes), 
but the residuals exhibit what appears to be a parabola. The 
residuals, however, seem to show a periodic excursion on either 
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side of the parabloic fit, with amplitude about 20 minutes and 
period about 35 years. Note that the residuals can also be fit with 
a series of five straight lines. These are the two paradigms men- 
tioned in the introduction. Note too, that after HJD 53000, the 
residuals fall further and further below the parabola which fit so 
well for 70 years, but the observations continue along the fifth of 
the straight line segments. 

Fitting the timings with a parabola requires some care be- 
cause of the periodic excursion on either side of that curve. It is 
thus important to use only a whole number of cycles to minimize 
the effect of the excursions on the parabola. The years 1934 to 
2001 meet this requirement and have adequate timings; when a 
quadratic least squares fit is made to just those years, the result- 
ing ephemeris is 



HJD miX = 28802.5661 



• 0.20102934 E + 0.755 x 10" 10 E 2 



(2) 



which is HGF Eq. 6 (note that their equation has an error in 
the quadratic term exponent). Fig. [2] shows the residuals after 
removal of the parabolic fit from Eq. 2. It also shows the fit to 
those residuals presented in HGF Eq. 8: 



(O-Qum = A [ 



0.7791 sin(<D+5.274) 
1.0+0.47 cosO 



- 0.3986 ] 



(3) 



where A = 0.013630223 in days and <D = 0.000 10084(£ - 7934) 
in radians. This function is plotted in Fig. [2] and when sub- 
tracted from those residuals, the new ones are plotted in Fig. [3] 
In both of those figures, the increasingly poor agreement after 
HJD 53000 can easily be seen, and amounts to more than two 
hours by 2012. 
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Fig. 2. Residuals for HGF Eq 6 (my Eq 2), a quadratic fit to all 
data. The solid line is the residuals expected from LiTE, HGF 
Eq 8 (my Eq 3) 
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Fig. 3. Residuals after removing HGF Eq 8 (my Eq 3), the LiTE 
fit. 



If it is included, a period change will be masked for several ob- 
serving seasons because the least squares fit will automatically 
adjust to the new, slightly disparate timings. 

Fig. |4] shows the residuals of this ephemeris to the fitted 
points (diamonds) and the new points (squares) after HJD 52000 
(year 2000). Though the new points seem to have slightly pos- 
itive residuals, there is no net trend away from zero, so no new 
period change is indicated. While the continuation of the lin- 
ear ephemeris does not guarantee this second interpretation to 
be correct, it is certainly consistent with it. This brings up the 
question of how the transition in period takes place. 
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Fig. 4. Residuals from the linear fit given in Eq 4. 



The other paradigm for the period variations in BW Vul, pro- 
posed by Tunca dl978l l and amplified by Chapellier & Garrido 
( 11990b . is that the period remains constant for one or more 
decade, and then undergoes an abrupt period change. To test the 
new timings against this model, a linear ephemeris was fit to the 
data from 1982 to 1998. which yielded an ephemeris: 

HJD^ = 28801.9903 + 0.2010439 E (4) 

Recall that the cycle count for this equation is still based on a 
cycle count where the Huffer 1938 maximum is considered to be 
cycle zero. It is important, when applying an ephemeris function 
to new data to not include the new data in the least squares fit. 



Chapellier & Garrido (1 19901 1 analyze the 1981 period change 
and conclude that it showed a transient phase over the years 
1979 to 1982 (about three years), with cycle-to-cycle variations 
of both amplitude and phase. However, they present only four 
sparse light curves, one consisting of only six measurements. 
Their Fig. 2 shows a large scatter in the residuals from 1981, but 
this was based on just one measurement. 

Fortuitous good luck brought to light additional data ob- 
tained by James Kemp in 1981 when he was doing polarimetry 
of BW Vul, but neither the polarimetry nor the concommitant 
photometry were ever published. The photometry is not optimal, 
as a rather faint nearby star was used as a comparison, but it does 
allow eight new timings (listed in Table [2]) to be derived for a 
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Table 2. Timings for BW Vulpeculae from James Kemp 1981 



Cycle 


HJD (2400000+) 


Max/Min 


(o-cy 


79793. 


44843.8917 


Max 


9.9 


79798. 


44844.8982 


Max 


11.8 


79803. 


44845.8962 


Max 


1.4 


79807. 


44846.7097 


Max 


14.8 


79813. 


44847.9131 


Max 


10.7 


79818. 


44848.9161 


Max 


7.5 


79887. 


44862.7868 


Max 


5.6 


79887.443 


44862.8838 


Min 


17.0 



(O-C) 1 is the residual to the linear ephemeris of Eq. 4. 



an inordinate fraction of (3 Cep stars in the shell-burning phase, 
as pointed out by Jerzykiewicz ( 11999b . 

It is anticipated that another period change should take place 
in the not-too-distant future, so continued monitoring of the tim- 
ing and light curve could prove to be very important. 
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critical year in the transition that previously had only one timing. 
These observations were communicated by Dr. Gary Henson. 

Fig. [5] shows an enlarged and extended view of Fig. ^be- 
tween 1976 and 1984. The diamonds are part of the timings that 
were fit by Eq. 4; the squares were not fit at all, but lie along a 
straight line which represents the shorter period prior to 1980; 
and the triangles represent timings in the 1980, 1981, and 1982 
seasons (also not fit). While it seems as though they lie on a curve 
above the straight line fits, they are only above by a few minutes, 
and it isn't clear if this is significant; certainly the scatter is al- 
most as large as the offset itself. Further, the scatter in those years 
is comparable to the scatter in prior and subsequent years. There 
doesn't seem to be the instabilities in period that Chapellier & 
Garrido (1990) have claimed. Their suggestion that the ampli- 
tude was also varying by 15% could be due to inhomogeneous 
data from different filters as well. There was no variation at that 
level seen in the light curves of the international campaign of 
Sterken et al (1987), which was also in the transition years. 
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Fig. 5. Residuals from the time of the 198 1 period change, based 
on the ephemeris in Eq. 4 
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4. Conclusions 

Photometry of BW Vulpeculae since 1995 fails to verify the 
continued secular period change or the cyclic period variation 
attributed to the LiTE of a binary system. It appears that the 
proper interpretation of the period variation in this star is a piece- 
wise linear ephemeris, ie constant period interrupted each one 
to three decades by an abrupt period change. If the discrete pe- 
riod changes are not evolutionary, this allows the star to be in 
the core hydrogen burning phase of evolution, consistent with 
the position of cluster members of this class of star in the color- 
magnitude diagram. It thus seems to help solve the problem of 



4 



